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Abstract: We present a detailed statistical study of the observed anisotropy in radio 
polarizations from distant extragalactic objects. This anisotropy was earher found by Birch 
(1982) and reconfirmed by Jain &: Ralston (1999) in a larger data set. A very strong signal 
was seen after imposing the cut |RM — RM| > 6 rad/m^, where RM is the rotation measure 
and RM its mean value. In this paper we show that there are several indications that 
this anisotropy cannot be attributed to bias in the data. We also find that a generalized 
statistic shows a very strong signal in the entire data without imposing the RM dependent 
cut. Finally we argue that an anisotropic background pseudoscalar field can explain the 
observations. 
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1. Introduction 



Polarizations of radio waves from extragalactic sources undergo Faraday Rotation upon 
propagation through galactic magnetic fields. This effect provides very useful information 
about astrophysical magnetic fields (Zeldovich, Ruzmaikin &; Sokoloff 1983, Vallee 1997). 
The amount of rotation is proportional to the magnetic field component parallel to the 
direction of propagation of the wave and to the square of the wavelength A. The observed 
orientation 0(A) of the linearly polarized component of the electromagnetic wave can there- 
fore be written as, 

e(A2)=x+(RM)A2 (1) 

where the slope, called Faraday Rotation Measure (RM), depends linearly on the line inte- 
gral of the parallel component of the magnetic field along the direction of propagation of 
the wave and x is the intercept, also called the intrinsic position angle of polarization, IPA. 

The observed polarization angle Xi after the effect of Faraday rotation is taken out of 
the data, is observed to be dominantly aligned perpendicular to the orientation axis of the 
galaxy. Let ^ denote the orientation angle of the galaxy. Then the offset angle |3 = x ~ 4" 
is found to be approximately equal to n/2 for most of the sources, i.e. the distribution of (3 
over a large sample of sources is found to peak at (3 ~ n/2. Besides this dominant trend, a 
smaller peak is also found at |3 « 0. This suggests the existence of two populations, some 
with polarization position angles parallel and others perpendicular to the galaxy axis. 

In 1982, Birch empirically observed an angular anisotropy in the offset angle (3, using 
a data set of 137 points. Birch's statistics were questioned by Phinney & Webster (1983) 
and it was pointed out that the significance of Birch's result can be significantly reduced 
if the experimental errors in |3 are taken into consideration. Phinney & Webster (1983) 
also suggested that the signal observed by Birch might result from the presence of bias in 
data. Kendall & Young (1984) further investigated Birch's claim of cosmic anisotropy with 
more sophisticated statistics and using an updated version of Birch's data. They found that 
the statistics were not consistent with isotropy at 99.9% confidence level. Later Bietenholz 
&; Kronberg (1984) repeated the calculations using single- number correlation test statistic, 
originally proposed by Jupp & Mardia (1980). This statistic also showed strong evidence 
of anisotropy in Birch's data with a confidence level of 99.9S%. They went on to create an 
independent set of 277 points which, however, showed no signal of anisotropy. This lead to a 
dismissal of Birch's results but left unresolved the puzzling fact that his data had contained 
a signal of anisotropy at such a high level of statistical significance. 

The possible existence of anisotropy in radio polarizations was reanalysed by Jain & 
Ralston (1999). The authors collected an independent set of 361 points from all the available 
catalogues. This set included (3 values for 29 sources which were contained only in the 
Birch's compilation. Jain & Ralston (1999) also considered the data set of 332 sources 
obtained after deleting these 29 objects. The authors found that both the data sets showed 
a statistically significant signal of anisotropy. The observed signal can be expressed as 
follows, 

2|3-7Tf«A-R (2) 

where R is a unit vector in the direction of the source. The vector A represents the three 
parameters of this fit and points in the prefered direction of this dipole anisotropy. The 
unit vector A controls the distribution of (2(3 — tt) on the dome of the sky such that it is 
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predominantly positive in the direction of the axis while it is predominantly negative in the 
opposite direction. 

Jain & Ralston (1999) tested the isotropy as null hypothesis using Maximum Likelihood 
analysis by taking into account the transformation property of (3 . A typical null distribution 
for angular variables is the von Mises (vM) distribution 

vM(e) = const X e'^=°'^(^' (3) 

where k is a constant and 9 is the angular parameter which is taken here to be 2(3. The 
factor 2 arises since (3 refers to the polarization angle and |3 + tt is identified with (3 (Ralston 
& Jain 1999). The joint distribution of the angle (3 and source positions was taken to be of 
the form 

f(2|3)C(2|3,R)h(R) , 

where f(2(3) is the marginal distribution of (3 which is assumed to be the von Mises, h(R) 
is the angular distribution of the sources and the correlation ansatz is taken to be of the 
form, 

C(2|3,R) =exp[A-Rsin2|3] (4) 

Here A represent the three parameters of the correlation ansatz. An alternate distribution, 
obtained by replacing 

2|3 ^ a((3) =2|3 +-vsirL2|3 (5) 

was also considered by Jain & Ralston (1999) since it was found to provide a better null fit 
to the data. The fit selected a negative value of the parameter -v ~ — 1 which leads to a more 
sharply peaked distribution in comparison to the von Mises. Such a sharply peaked function 
is also justified by the distribution of polarizations obtained under the assumption that the 
electric field components of the electromagnetic wave are normally distributed (Sarala & 
Jain 2001, 2002). 

Likelihood analysis showed a strong signal of anisotropy using the transformation [S] in 
the joint distribution \^ The strongest signal was obtained after making the cut on the 
rotation measures so as to eliminate the data with RM in the range |RM — RM| < 6, where 
RM = 6 rad/m^ is the mean value of the rotation measure over the data sample. In fig. 
^we show the histogram of the rotation measures for the 332 sources. It is found that the 
peak position as well as the mean value of RM is shifted from zero and lies approximately 
at 6 rad/m^. This is somewhat unexpected since for a large unbiased sample we would 
expect the peak value of RM to be zero. Even for a larger data sample compiled in Broten, 
MacLeod & Vallee (1988) we find similar distribution with peak position as well as the mean 
and median of the rotation measures shifted from zero towards positive values. The cut 
on RM, therefore, removes points lying close to the central peak in the distribution. The 
statistical significance of the signal after this cut was found to be approximately 0.06%, i.e. 
the probability that the correlation seen in data might arise in a random sample is given by 
P = 0.06%. The unit vector A after making the cut |RM — RM| > 6 was found to be (Jain 
& Ralston 1999) 

^ = [(0 h,9 m) ± (1 h,Om),-r± 15°] (6) 

in the B1950 coordinate system. 

Similar results were obtained by using the Jupp and Mardia correlation coefficient. In 
this case the authors (Jain & Ralston 1999) found the correlation between sin(2(3 — sin2(3] 
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and the angular position vector of the source. The parameter "v in Eq. [21 was set equal 
to —1 since the likelihood analysis generally prefered a value between —0.5 and —1. The 
statistical significance in this case, after making the cut on rotation measures, was found 
to be P = 0.04%. 

In the present paper we further illustrate the nature of this correlation. The basic cause 
of this correlation is so far unknown. There have been several proposals such as observa- 
tional bias, bias in the extraction of RM and (3 (Phinney & Webster 1983) existence of a 
light pseudoscalar (Jain, Panda & Sarala 2002), cosmic rotation (Obukhov 2000, Obukhov, 
Korotky & Hehl 1997, Mansouri & Nozari 1997), inhomogeneous and anisotropic universe 
(Moffat 1997), weak lensing by large-scale density inhomogeneities (Surpi &: Harari 1999), 
background torsion field (Dobado & Maroto 1997) etc. The fact that the correlation is 
improved significantly after making a cut based on rotation measure may provide a hint 
about its origin. One obvious possibility is that there may exist bias in the extraction of RM 
and |3 which leads to a correlation between these two parameters. Then, since the RM or 
the magnetic field is not isotropically distributed, it can induce an anisotropy in p. Indeed 
such a possibility was suggested (Phinney & Webster 1983), soon after the appearance of 
the Birch's paper. 

2. Statistical Analysis of the Angular Correlation in Polariza- 
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Figure 1: The distribution of Rotation Measure for the data with 332 points excluding some 
of the sources listed in Birch (1982) for which information on RM was not available. The 
RM distribution is not exactly Gaussian. Its mean as well as the peak position is shifted 
from RM = towards positive values, as discussed in text. The vertical dashed line shows 
the position of the mean, RM = 6. 
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tions 



We next perform several statistical tests in order to determine whether the observed 
correlation can be explained in terms of bias in the values of (3 which may result from 
biased extraction of rotation measures. As discussed above if the extraction of RM and (3 is 
biased then since the RM is not isotropically distributed on the celestial sphere it can lead 
to an anisotropic distribution of (3. In order for this mechanism to work it is necessary that 
(3 is strongly correlated with RM. We next examine if such a correlation exists in the data 
by computing the Jupp & Mardia (1980) (JM) correlation between RM and (3 for the data 
set containing 332 sources. The correlation is found to be very small with np^ = 1 .742. We 
tried several cuts to see whether the correlation is enhanced and found a small correlation 
with np^ = 4.341 for |RM| < 100 rad/m'^ with the number of sources n = 278. A more 
stringent cut on RM does not improve the correlation. We summarize JM correlation with 
different cuts in Tabled We see from this table that RM is not correlated with (3. This 
implies that a large range of explanations of the anisotropy which invoke the possibility 
of biased estimation of RM and IPA, including the mechanism proposed by Phinney & 
Webster (1983), are disfavored. 
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Figure 2: Aithoff-Hammer equal-area plot of sin 2(3, where |3 is the difference in the polar- 
ization position angle and the orientation angle of the source. We represent sin 2(3 < 0.0 
by blue pluses while sin 2(3 > 0.0 is represented by red crosses. The dipole anisotropy axis, 
[(0 h.,9m)±(l h,0 m),-1°±15°], obtained by Jain & Ralston (1999) is shown along with 
its error bars. The position of the Virgo supercluster is shown with a red triangle. The cut 
RM — 6| > 6 has been used in this plot. 

The above lack of correlation between RM and (3, however, does not rule out the possi- 
bility that the value of (3 is randomized due to biased extraction of RM. This might happen 
preferentially in the regions of large magnetic field i.e. where RM is large. As we have 
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discussed earlier, the distribution of (3 peaks at 7t/2. If in a certain patch the values of (3 
are randomized, then the mean value of (3 in this patch will deviate significantly from 7t/2 
and could in principle give rise to a large value of the statistic indicating correlation. 

We next investigate whether the correlation arises due such a randomization. In Figure 
121 we plot the sign of sin2|3 over the dome of the sky in the equatorial coordinate system 
after making the cut |RM — 6| > 6. Here the sources with sin2|3 < 0.0 and sin 2(3 > 0.0 are 
represented by plusses and crosses respectively. From this figure we see that the positive 
values of sin 2(3 dominate near the dipole anisotropy axis given in Equation |BJ Figure [3 
clearly indicates that there indeed exists an anisotropy in the polarizations of radio galaxies 
and there exist large patches where the value of sin 2(3 is either positive or negative. Hence 
we find no evidence that the correlation seen in data arises due to randomization of (3 in 
some regions. We may quantize this clustering of negative and positive values of sin 2(3 
by defining a new statistic. We consider a Nn number of nearest neighbours of the source 
i. For these sources we set X| = 1 if sin 2(3 > and X| = —1 if sin 2(3 < 0. We 
then compute Xt = | ^^xtl and the statistic X is obtained by summing X^ over the entire 
data sample. The P-value in this case is computed by comparing with a large number of 
random samples obtained by shuffling the (3 values of the data. The resulting P-values or 
the significance level for the complete data sample and for the set obtained after making 
the cut |RM — 6| > 6 are shown in Figure 01 We find that P ~ 0.2 — 0.5% for the complete 
sample and P ~ 2.5 x 10^ — 6.0 x 10^"^% after the cut. In Figure^lwe show the histogram 
of the statistic X generated using 10000 random data samples. This figure also clearly 
shows that the statistic of the data sample lies far above the peak value of the histogram. 
The very small P-values obtained after the cut is a clear indication that (3 values in the 
different regions of the sky are indeed correlated with one another. This is infact another 
confirmation of the large correlation found by Jain &: Ralston (1999). Hence we dismiss the 
above explanation in terms of bias in the extraction of IPA. 

An improved statistical procedure for extraction of rotation measures and IPA was 
proposed by Sarala & Jain (2001, 2002). Using the revised (3 values obtained by this 
procedure, the Jupp-Mardia test for anisotropy gives P = 0.1 %. Hence these revised (3 
values continue to support the presence of anisotropy in the data. 



Cuts 


n 




Full data 


332 


1.742 


IRM-61 > 6 


265 


1.780 


|RM| < 100 


278 


4.341 


|RM| < 70 


254 


1.101 


IRMj < 100 and |RM - 6| > 6 


211 


4.320 



Table 1: JM correlation between RM and (3 with different cuts in the RM. 

In order to further examine the effect of rotation measures on the observed anisotropy 
we examine how the correlation changes if we eliminate the sources which lie close to the 
galactic plane. The dominant contribution to RM is obtained from the milky way and hence 
largest values of RM are observed for these sources. We therefore impose a cut on the data 
to eliminate sources with the galactic latitude |b| < 30°. This leaves a total of 214 sources in 
the data set. We repeat the JM correlation discussed earlier using Map 3. We find that the 
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Figure 3: The P-value(%) using the X statistic is plotted as the function of the number of 
nearest neighbours. The result is shown for both the complete set(332 points) and with the 
cut in the RM (265 points). 

statistic np^ = 9.43 after the galactic cut and this imphes that the anisotropy in (3 persists 
with about 98% confidence level. As expected, the anisotropy is enhanced by making a 
cut in the RM such that only the sources which statisfy |RM — 6| > 6.0 are included. The 
results are tabulated in Table l2j We see from the table that the statistic does not change 
much in comparison with the data without the galactic cut. Hence we find that the regions 
of large RM do not necessarily imply large correlation of (3. 





n 




P(%) 


Full Set 


214 


9.43 


2 


|(RM-6| > 6.0 


155 


22.33 


0.04* 



Table 2: The np^ values after eliminating sources lying in the galactic plane. Results are 
given for the full set as well as with a cut in the Rotation Measure. *The P value for the 
cut |(RM — 6)1 > 6.0 is obtained by using Xs- 

The above discussion clearly shows that the correlation seen in the data is not a direct 
consequence of the correlation of the RM with the angular coordinates of the source. How- 
ever the cut in RM does result in a large increase in the signal of anisotropy and hence in 
some sense it must depend on RM. One possibility is that the sources eliminated by the 
cut show a random behaviour and are uncorrelated. This turns out to be not true. Instead 
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Figure 4: The histogram of X statistic using 10000 random samples obtained by shuffling 
the (3 values of the data with the cut |RM — 6| > 6. The number of nearest neighbours is 
taken to be 135 at which the P- value is found to be the smallest. The statistic of the data in 
this case is equal to 0.165, shown by the dotted line, and the corresponding P = 2.5 x 10^ 
%. 

we find that the sources eliminated by the cut show an opposite angular dependence in 
comparison to the remaining sources. We illustrate this by making the following ansatz for 
the joint distribution 

F(|3,R)li(R) =exp[kcosa + A-R(l + g(RM)) sinO] li(R) (7) 

where as defined earlier O = 2(3 +'vsin2(3 and g(x) is taken to be a gaussian, 

g(x)=ci7Texp -[x-czf/cl . (8) 

Here we have introduced three new parameters Ci,C2 and C3. The difference of maximum 
likelihood between this correlated ansatz and the null ansatz turns out to be 17.6 for the 
entire data set of 332 sources. Here the null ansatz is the two parameter distribution 
obtained by using the transformation[21in the von Mises distribution. This "sharply peaked" 
distribution can be written as 

S(|3) = const X e'^=°''(2|3+-vsin2|3) (g) 

The statistical significance of the signal of correlation in this case is 

P = 3.8 X 10^4% , 
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taking into account the eight parameters in the correlation ansatz. If the null hypothesis 
is assumed to be the von Mises distribution we find an even smaller value of P. The signal 
of correlation is clearly very strong and cannot be dismissed as a statistical fluctuation. 
The parameters are given by k = -0.561 ± 0.071, v = -1.09 ± 0.17, Ci = -1.71 ± 0.32, 
C2 = 9.64 it 0.47, C3 = 2.89 it 0.55, |A| = 0.70 ± 0.15 and the axis parameters are given by 

RA = (0^, 3.7"^) ± (0^, 46""), DEC = -10.7° ± 12.2° (10) 

in the B1950 coordinate system. Within errors this points opposite to the direction of the 
Virgo supercluster. We see from the value of the parameter c-[ that the sources whose RM 
lies roughly within the interval C2± C3 show a correlation with A • R which is opposite to the 
remaining sources. The correlation term A • R(1 + g(RM)) has an opposite sign for these 
sources in comparison to the remaining sources. The axis parameters, Eq. 1101 obtained 
from this fit agree with the earlier fit, Eq. El within errors. 

The function g(RM) essentially treats the sources whose RM lies in the interval C2± C3 
differently from the remaining sources. This interval is centered at the value Cj which is 
very close to the position of the peak in the RM distribution as can be seen in Fig. ^ 
We may attribute this shift in the peak to some bias in the extraction of RM or due to 
some systematic physical effect such as the contribution due to the milky way. Hence if the 
difference RM — 02 is small it may indicate that the host contribution to RM is negligible. 
The introduction of the function g(RM) in the correlation ansatz, therefore, treats the 
sources with negligible contribution due to the host galaxy differently than the remaining 
sources. 

The relationship between |3 and the position vector R of the source indicated by Eq. [3 
may be expressed as follows, 

<n+.cos2p)sinQ> ^^_R 

< (1 +-vcos2(3jcosO > k 

where the angular brackets indicate statistical averages. Since the parameter v ~ — 1, the 
function (1 +'y cos 2(3 ) is approximately 2 near the peak of the distribution of (3 and zero near 
the tail. Hence this function weights the peak region more than the tail in the statistical 
averages. For |3 close to its peak value n/2, we find that sinO ~ —2(1 — ^)(|3 — 7t/2) and 
cos O ~ — 1 . If RM is not too close to its peak value or more precisely |RM — C2I > C3 then 
g(RM) ~ 1 and hence we find the approximate relation 

» 1 A • R , , 

which is same as the relationship given in Eq. [21 with the A in Eq. (21 identified with ^^^^^^ 
in Eq. HH 

The anisotropic distribution of the offset angle obtained with the fit Eq. |3 may now 
be illustrated by making an Aithoff-Hammer scatter plot of the correlation variable (1 + 
g(RM))sinO. This is shown in fig. where stars and plusses represent sources for which 
(1 + g(RM))sinn > 0.1 and (1 + g(RM))sinn < 0.1 respectively. We find that the stars 
are concentrated in the direction of the axis with the plusses concentrated in the opposite 
direction. The anisotropy in the angular distribution of the offset angles is clear from 
this plot. The nature of the anisotropy is also nicely illustrated by making a scatter plot 
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Figure 5: Aithoff-Hammer scatter plot of the correlation variable (1 + g(RM))sinQ, 
as defined in Eq. [7| Here red stars and blue plusses represent sources for which 
(1 + g(RM))sinO > 0.1 and (1 + g(RM))sinn < —0.1 respectively. The anisotropy axis, 
[(0'\ 3.70™) ± (O"^, 46™),— 10.7° it 12.2°], along with its error bars is shown in green. 

after taking the local averages. In fig. |21 we show the Aithoff-Hammer scatter plot of 
the averaged correlation variable, (1 + g(RM))sinQ, at the location of each source. The 
average is done over the nearest neighbours whose angular separation from the source is 
less than cos^^ 0.95 ~ 18°. The mean number of sources which contribute to this average is 
found to be 12.6. The plot clearly shows anisotropic nature of the distribution of the offset 
angles. 

3. Physical Origin of the Effect 

From the results obtained above we see that the offset angle |3 shows some correlation 
with RM but the correlation is direction dependent. It is difficult to visualize how such a 
direction dependent correlation can arise due to bias in data. In all likelihood this is caused 
by some physical phenomenon. The offset angles of the different sources may be intrinsically 
correlated or the effect may arise due to propagation. The possibility of intrinsic correlation 
can be further investigated by determining the correlation of the intrinsic position angle of 
polarization x ^-iid the orientation angle of the galaxy ^\) among the different sources. By 
using the statistical techniques used in (Hutsemekers 1998, Hutsemekers &: Lamy 2001, 
Jain, Narain & Sarala 2003) to test of alignment of optical polarizations, we find that 
X and \|j for different sources show no alignment with one another. The P value in this 
case is found to be larger than 10 % independent of the number of nearest neighbours 
used to test for alignment. Hence it is not possible to attribute the radio anisotropy to 
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intrinsic alignment of radio polarization and it is likely to arise due to some propagation 
effect. Within conventional physics the corrections to the Faraday polarization rotation 
effect are very small and hence it is unlikely that this anisotropy may be explained within 
this framework. The anisotropy may, therefore, be an indication of some new physical effect. 

One interesting possibility is the presence of a very low mass background pseudoscalar 
field in the universe. We assume that at large redshifts, relevant for the sources considered 
in the present data set, this field has a coherent component which is approximately given 
as 

ct)(R) «c|)ocos(e)(l +g[B,(R)]) (13) 

where (^o is a constant independent of redshift, 9 is the angular position of the source 
measured with respect to the prefered axis A, B,- is the radial component of the magnetic 
field in the coordinate frame with origin located at the observation point i.e. earth and glBr) 
is a gaussian defined in Eq. |HJ At the position of the earth we assume that this field acquires 
a value cj)!. The total rotation in the polarization angle due to this field is then equal to 
ct)(R) — (Harari & Sikivie, 1992). Since the rotation measure is proportional to B,- we 
find that this effect can explain the observed correlation given in Eq. [3 This mechanism 
requires a large scale anisotropy in the universe and provides the simplest explanation of 
the observations. 

We next determine whether the observations can also be explained within the frame- 
work of an isotropic universe. Although the correlation indicates the presence of a global 
anisotropy it is possible that the correlation may arise due to the existence of a few patches 
of length scales of order Gpc over which the pseudoscalar field shows a coherent dependence 
as a function of the angular coordinates. For example there may exist two such patches, one 
in the direction of the axis and another in the opposite direction. The observed anisotropy 
can arise if the pseudoscalar field in one of these patches is approximately equal to the 
negative of its value in the other patch. In particular the observed correlation can be ex- 
plained if ct)(R) ^ 4>o(l +g(l^M)) in one of the patches and c|)(R) ^ — ct)o(l +g(RM)) in the 
second patch. We point out that although the existence of a large scale dipole distribution 
of the background field violates the fundamental assumption of isotropy of the universe, the 
existence of the few patches of length scales of order Gpc does not. Here we also assume 
that in eq. El4)(R) depends on the magnitude |B| rather than on Br. Such a functional 
dependence may also explain the observations provided we attribute the RM contribution 
due to the host galaxy to the deviation of RM from its peak value, RMpcak, in the RM 
distribution. As discussed earlier, our fits also always select the parameter Cj such that it 
is approximately equal to RMpoak- The shift in RMpeak from zero may be attributed to 
some local effect such as the contribution due to the milky way or due to some bias in the 
extraction of RM. The deviation of RM from its peak value, RM-RMp^ak; gets contribution 
both from the host galaxy and the milky way. Hence RM-RMpeak is directly proportional to 
the magnitude of the magnetic field in the host galaxy, assuming that in a statistical sense 
all the components of the background magnetic field have equal strength. This mechanism 
may explain the observed correlation within the framework of an isotropic universe provided 
their exist patches of coherent background pseudoscalar field of length scales of order Gpc. 
However coherence on such large scales is not expected in conventional cosmology. Fur- 
thermore this explanation also requires an accidental correlation of two oppositely directed 
patches such that pseudoscalar field is positive in one of these patches and negative in the 
other. 
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A large scale anisotropy has also been observed in the optical polarizations from distant 
quasars (Hutsemekers 1998, Hutsemekers & Lamy 2001, Jain et al. 2003). It was found 
that the optical polarizations are aligned with one another over very large distances of the 
order of Gpc. A very strong alignment was seen in a patch at large redshifts centered at the 
Virgo supercluster (Hutsemekers 1998). Furthermore it was found (Jain et al 2003) that the 
polarizations of the large redshift, z > 1 , data sample are aligned over the entire celestial 
sphere. This effect is also not easily explained within conventional cosmology /astrophysics. 
The effect may have some relationship to the radio anisotropy discussed in this paper since 
both the effects appear to be very strong in the direction of the Virgo supercluster. As 
shown in Jain et al. (2003) the optical alignment effect can also be explained if we assume 
the existence of a light pseudoscalar. The proposed explanation requires the existence of 
several patches of large scale magnetic fields of length scales of order Gpc. Furthermore the 
global alignment at large redshifts requires the existence of a magnetic field coherent over 
the entire universe. Alternatively this very large scale alignment can be explained if the 
large redshift QSO's emit a significant flux of light pseudoscalars at optical frequencies. 

4. Conclusions 

To conclude, we find that there is considerable evidence for the presence of angular 
correlation in the radio offset angles which is not easily explained within conventional cos- 
mology/astrophysics. The effect may be explained by the presence of a hypothetical pseu- 
doscalar particle. The effect indicates the existence of a global anisotropy in the universe. 
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Figure 6: AithofF-Hammer scatter plot of the averaged correlation variable x = 
(1 +g(RM))sinO, as defined in Eq. d For each source located in the direction Rt we 
take the average over all the nearest neighbours which lie within the cone defined by 
Ri • R2 < 0.95, where R2 is the unit vector in the direction of any other source. The 
anisotropy axis, [(0*^, 3.70™) ± (0*^, 46™),— 10.7° it 12.2°], is shown along with its error 
bars. Within errors it points opposite to the Virgo supercluster. 
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